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1. Introduction. 

The Hopf map gives the structure of the three-sphere as a fibering, with a two 
sphere as base and a circle as fibre. The projection — > can be expressed in 
terms of the Cartesian coordinates on the corresponding embedding and by 
the non-hnear mapping (Hopf [1]), 

yl = 2{x^x^ ±x^x^) (1) 

The puUback of a harmonic polynomial in a set of ys is a harmonic polynomial 
in the xs. In this paper I wish to expand on this fact, which is useful because it 
relates eigenfunctions on and S^. These two facts are really the same. 

Kibler et al, [2], have used the Hopf mapping to relate modes. Hage Hassan 
and Kibler, [3] , have elaborated on this aspect using the Fock-Bargmann-Schwinger 
approach to angular momentum theory. Some of the material I give in the next 
section will be found in these works, and elsewhere. 

A fibre bundle treatment is given by Boiteaux, [4], (see also Gilkey et al, [5]) 
but my discussion will be more simple minded. 

In later sections I extend the discussion to include homogeneous factors of 
and will give a practical means of obtaining the modes. In this regard, the Hopf 
approach has also been expounded by Weeks, [6], and Lachieze-Rey and Weeks, [7], 
in an astrophysical context and I refer to these works for references to some other 
mode calculations. I should also mention the interesting work of Bellon, [8], [9]. 

An important step on the road to modes on the factored three-sphere is the 
construction of modes on the orbifolded two-sphere. These polyhedral harmonics, 
also known as crystal, or lattice, harmonics, have been investigated for many years; 
early papers being by Elert, [10], and Bethe, [11]. In particular, more recently, 
icosahedral harmonics have found use in the study of viruses, fuUerenes and quasi- 
crystals. The method presented here gives a relatively simple way of finding them 
and information in [8] allows the tensor icosahedral harmonics to be found too. 
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2. Harmonic polynomials and eigenmodes 

Working with polynomials has its attractions, but the puUback result is most 
easily appreciated in angular coordinates. The Euler angle expression for the a;*, 
going back at least to Klein and Sommerfeld, is (I repeat some standard relations),^ 

x'^ = -Rsmd/2 sm(ip - (/))/2 
x"^ = Rsme/2 cos(V' - 0)/2 

x^ ^Rcose/2sm{ilj + (f))/2 ''^^ 
= i?cos6'/2 cos{ip + (j))/2 , 

whence, from (1), the V'-independent, and (/)-independent, combinations, 

y]^ = r sin 9 cos (/> , y\ — r sin 9 cos ip 

y\ = r sin 9 sin (j) ^ y^ = r sin 9 sin ip (3) 
y^ = r cos 9 , y^ = r cos 9 

with r = R^, implying the map S*^ ^ S^. In fact, lifts to a quartic surface in M"^ 
which is the product of two (identical) three- spheres. 

I have chosen the factors so that the SU(2) element, U, agrees in form with 
that in [12], equn.(2.21), viz, 

U = = a-oj _ ^x.o- (4) 

in terms of Pauli matrices. 

Each set of ys in (3) produces a system of spherical polar coordinates on a 
corresponding R*^. I will use just the y_^, and denote these by y*. A complete set 
of Laplacian cigcnfuctions on is provided by the spherical harmonics, Cl^{9,(p). 
The relation between these and the SU(2) rep matrices, T>-l^^{9, 4>,ip)j is {e.g. Brink 
and Satchler, [13], Vilenkin, [14], Talman, [15]), 

cr(^,0) = (-i)^pi^°(^,0,v), lez, (5) 

the right-hand side being independent of the fibre angle, Since the Vs are 
Laplacian eigenfunctions on ~ SU(2) the result follows. It is standard {e.g. Hill, 

2 My Euler angles correspond to the ''z-y-z^ rotation convention, adopted by [12], [13]. Note that 
Vilenkin, [14], uses the other popular choice, z-x-z. 
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[16], [12]), that the sohd quantity, i?^'©', is a harmonic polynomial in the of 
degree 21 and the corresponding eigenvalues equal 4Z(Z + 1) agreeing with the 
eigenvalues from (5) taking into account the scaling relation between the relevant 
Laplacians (see below). 

This can be made more explicit by introducing the relevant differential opera- 
tors such as the right-generators of the SO (4) action on S^, 

y_i = Y,* (6) 

the first two of which allow the right index in (5) to be raised and lowered by the 
standard action, 

-iY^V^{g) = V^{g)J^, (7) 

where are the spin-j angular momentum matrices. 

The Laplacian, Ag, on the unit three-sphere is the Casimir operator. 



A3 = 4Y^ = 4{Y,Y_, + Y_,Y, - Y^^) 

= 4 ( a| + cosec ^e{dl + dl) -2cos9 cosec ^dd^d^ + cot 9 



(8) 



This corresponds to the metric,^ 

ds^ = ^ {de^ + d(p^ + dilj^ + 2 cos 9 d(p dip) 
= ^ {d9'^ + sin^ 9 d(j)'^ + (dV + cos 9 d(j)f) , 
which occurs as the spatial part of the metric on R^, 

df = dR^ + i?2 ds^ , 
giving the Dalembertian (see e.g. Hund, [18]), 

□4 = -^3 (i?^ a^) + A3 . 



^ The numerical factor can be checked by computing the volume of the three-sphere using the 
conventional ranges O<0<7r,O< 0<27r,— 27r<V'<27r (but see [12], [17]), and the metric 
determinant ^ = sin 0/8. The volume is 27r^o^ for radius a. cf [5]. 
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Now let ^4 act on a puUback, i.e. on a ■^-independent function, fQ(9,(f)), e.g. 
Cp(9, 0). Then from (8), and using r = R'^ and Yq /q = 0, 

= 4r(la.(r^a.) + lA,)/o 
= 4rn3/o, 

where A2 is the unit two-sphere Laplacian, in with coordinates y"', 

A2 = + cot ed0 + cosec dl . 

Equation (9) shows that the Hopf mapping pulls harmonic functions back to 
harmonic functions. Most of the above is well known, in one form or another. 



3. Harmonic projection 

It is relevant to consider, briefly, the relation between harmonic projections in 
the light of the Hopf map. 

A standard result is that, in an R^, a rational integral polynomial of degree n, 
/„(y), can be decomposed as 

/„(y) = yjy)+yVn-2(?)- (10) 

Y^if) is a harmonic polynomial (solid spherical harmonic). For d = 3 this was 
proved by Gauss and his method easily extends to (i-dimensions. 

The polynomial Y^ is the harmonic projection of and there is an explicit 
expression for it, {e.g. Hobson, [19], Vilenkin, [14]), 



1 - T T + 



2(2n + d-4) 2.4(2n + d-4)(2n + d-6) 



obtained from the iteration of (10) and applications of (c/ Clebsch, [20], for 
d = 2>). 

I apply this formula for d = ?> and d = A. In the former case, to accord with 
my previous notation, I use as Cartesian coordinates. I also write y.y = and, 
for d = A, set f- = B?. 
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For d = 3 

and 

H{fn) = 

while for d = 4 



fn{y) = H{Uy)) + r'fn-2{y), 



+ 



and 



2(2n-l) 2.4(2n-l)(2n-3) 



fn 1 



1 - 



(12) 
(13) 
(14) 
(15) 



2.2n 2.4.2n(2n-2) 

In the special case of the Hopf mapping, harmonic projection commutes with 
hfting. If hfts to ^2^, then H{f^) Hfts to H{F2^) and, if /^_2 Hfts to -F2n-4' 
one has -F2n-2 ~ -^^-^2n-4 t>y comparing (12) and (14). Therefore, in this instance, 
(14) becomes 

Fjj:) = iy(Fjj:))+i?4F,_4(3:), (16) 

and also one sees that (13) lifts to (15). A direct proof, using the relation between 
Laplacians, (9), is not obvious. 

The extension to higher dimensions, d, of Thomson and Tait's approach to 
spherical harmonics is straightforward on noting that l/y*^"^ is harmonic, ^ 0, 
with j:^ = (x^)^ + (x^)^ + . . . + (x'^)^. This is an ancient fact,^ and another expression 
for the harmonic projection is, [19], 



H{fn) = (-1)' 



.2n+d-2 



{d-2)d(d + 2)...{d + 2n-Ay 



/n(V)^ 



(17) 



In particular, in three dimensions. Maxwell's multipole expression of a general 
solid harmonic is, 



n 

n.(y,P(fc)) = c?'"+'n(Pw^)-' 



fe=i 



which depends on the n 3-vectors, p^/,), k = 1, . . . ,n. 



To include the point r = formally, it would be best to adopt a distributional, Green function 
approach. See e.g. Rowe, [21]. 
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4. Spherical factors 



As described in section 2, lifting the modes from produces only the integer 
spin modes on S^. A complete set would also include the half-odd integer Vs. These 
can be eliminated by dividing by a Z2 antipodal action (to give the projective 
three-sphere) and requiring periodicity. (Relatedly, the Hopf map is unchanged 
under parity, y — > —p.) Alternatively, the integer spin modes could be eliminated 
by choosing anti-periodicity, as allowed by the topology, [22] , [23] , [24] . 

The total set of integer-spin modes, V^^' , is obtained, as mentioned, by acting 
with the raising and lowering right operators on the puUbacks of the modes, (5), 
in familiar angular momentum fashion. Weeks, [6] , see also [7] , describes this process 
using the polynomial approach and refers to the index m' as the 'twist'. A more 
conventional name is weight, coming from invariant theory, via Lie-group theory. 
It is measured by the vertical operator, Yq. 

Since the answer is already known, there is no especial calculational merit 
in finding the modes on 5^/1^2 in this way but it does give them an SO (3), 
geometrical character, arising originally from the isomorphisms S0(4)~ SU(2)x 



It is a particular example of the relation between the spectral problems on 
the free action S"^/r' and the orbifolded S^/F where V is the double of F, noting 
that F' = Z2 when F = 1. The eigenvalue aspects of this relation have already 
been expounded and used in [25] and now the geometrical underpinnings are more 



It is only possible to find the modes in this puUback way for a general (say left) 
symmetry action, S"^/F', if F' contains an antipodal Z2. This applies to even lens 
spaces, and when F' = O' , T' , Y' . 

The modes on S^/F' can be obtained by symmetry adaptation. I define, to 
begin with, the left group average, or projection. 



SU(2)/Z2 and SO(3)-SU(2)/Z2 . 



apparent, c/ [6], [7], [9]. 




7'Gr' 



(18) 



= 2C Y.^L^^lWn\Q.<^.i>) 



7e r 
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The last line is now the preliminary symmetry adaptation of the modes on the Hopf 
0)-sphere base.^ I do not pursue this approach to compute these modes but will 
present a better method in the next section. 

In the derivation of the relation (18) between the two projections, I have used 
the fact that T = F'/Zg where Z2 is an antipodal action, the non-trivial element of 
which corresponds to a rotation through 2n and is the 'central' element, introduced 
by Bethe to generate the binary group from the pure rotation one. For I integral, 
it is equivalent to the identity. Hence the factor of two. 

As before, the complete set of (preliminary) modes on S"^/F' is obtained by 
acting on the right with the raising and lowering operators. For a particular eigen- 
value (depending on only the label /), the degeneracy is the product of the right 
degeneracy, which is the range of the right index, i.e. (2/ + 1), and the left degen- 
eracy, d{l; F), evaluated by cutting down the overcomplete set of the 4)1^^ modes to 
the minimal number, usually by diagonalisation. The results were given in [23] and 
derived using the results of Polya and Meyer, [26], who employed Molien's theorem 
rather than constructing the group average directly. 

Most easily, the (left) degeneracies can be obtained from those on cyclic groups 
by making use of the cyclic decomposition of a spectral quantity S, 

I I ^ (J q ^ 

where the group F has axes of order q. For those groups satisfying the orbit- 
stabiliser relation, |F| = 2Q'nq, 

Si^V) = \{y^S{J.^) - SiJ^S) ■ (20) 

Relation (19) is given by Meyer, [26], and proved in [27]. It is the simplest 
way of deriving the degeneracies, and also their generating functions, which satisfy 
(20) and are given in [26] (see also Laporte, [28]). For example, the cyclic (Z^) and 
octahedral (O) groups, the generating functions, defined by 

00 

^(a;F)^^d(/;F)a', 
/=o 

^ I make V act on the left so that the coordinates on the sphere are the conventional d and 4>. 
The coordinate ranges and boundary identifications are detailed by Jonker and De Vries, [17], 
and exhibit the twisted product structure of S^. 
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equal 

gia; ZJ = — — ^1^^ , gia; O) = ^ . (21) 

' 1' i-(jl-ai' ' ^ (1-(j4)(1-ct6) ^ ^ 

The construction of the modes themselves is not so straightforward. 

5. Scalar polyhedral modes 

The easiest part of determining the independent set of symmetry adapted 
modes on S ^r' is the right raising and lowering, which is routine. The difficulty lies 
in the construction of the symmetry adapted spherical harmonics on S^/F, a topic 
of great interest to physicists, chemists and biologists over a long period. A lot of 
expository and technical material exists which it is impossible to survey. There are 
various approaches. In some, the group average is performed element by element. 
In others, it is bypassed, or camouflaged. 

I will spend a little time developing a method of calculating the symmetry 
adapted modes after which the puUback modes on S^/V can be considered to be 
known by the process outlined earlier. 

Perhaps the method that is most readily automated is one described first by 
Hodgkinson, [29], and based on the expression (13). The same technique was de- 
veloped later, independently by e.g. [30], [31]. 

Hodgkinson presented a complete process, if somewhat sketchily,^ and com- 
plained at the 'very heavy labour' involved in the harmonic projection last step, 
(13). Nowadays this can be alleviated using symbolic manipulation, and the al- 
gorithm has been detailed, independently, with this in mind by Ronveaux and 
Saint-Aubin, [32]. 

The method depends on the existence of an invariant polynomial integrity basis 
whose evaluation is classic, going back, in this setting, to Klein and with much 
subsequent work. 

The process yields a basis for invariant harmonic polynomials in the Cartesian 
coordinates of an embedding R . In the present notation, these coordinates are the 
y-^, y^, of (3). These polynomials can then be converted to polynomials on S^/V 
via the relation (1) and the right raising and lowering operators applied to complete 
the computation of a full basis, if it is desired to be so explicit. However, the results 
would not be illuminating. For large angular momentum, Cartesian coordinates 

^ Hodgkinson considers only those modes even under reflections in the polyhedral symmetry 
planes. He also does not give any specific computations. 
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become unwieldy. Better for the purpose if the harmonic basis on S^/F could be 
expanded in (standard) spherical harmonics, (for fixed I), for then, using the 
relation (5) with the irrep matrices, the right raising and lowering is immediate and 
consists simply of replacing the right index by m', running from / to —I. This is 
because left and right are completely independent for homogeneous factorings. 

I will now elaborate on the computational method outlined by Hodgkinson 
based on the form (17) rather than on (13). 

The tesseral harmonics can be defined in the Thomson and Tait, [33], way, {e.g. 
MaxweU, [34], Hobson, [35], Honl, [36]), 



yr" = pr (cos e) e^"^"^ = n,^ r'+^ (-1)"^ — —, — - 

y-m ^ pm/ g^ ^-im^ ^ m ^1+1 1 

for I > m > 0, — y.y, and 

iv.. = (-.)^^^^)"^ 



(22) 



{l-m)\ 

Instead of the Cartesian components, I use spherical ones, defined by 

- V +iy^ .3 

y.---^, 71". yo = -'y' 

for the 'standard' components, and 

y^ = y-i , y~^ = 2/1 , 2/° = -2/0 > (23) 

for the contrastandard ones. Also yl = y_^ and j/q = — j/q. 

If the Cartesian polynomials are real then one will need the combinations, 
(unnormalised spherical harmonics). 



1 al-m / am am \ 1 

y,r = Pricos e) COS m(P = -iV. r'+^ ^ + (-1)"^^^ - 

^ ^ ' V 2 dy^Q-'^\dy'^^ ' dy'^^jr 

1 al-m / am am \ i 

y,T = Pr(cos 9) sinm0 = -N^^ r^+^ - {-l)"^^— - 

^ ^ 2z dyl-"'\dyY' dy^H-^^J r 



(24) 



These are the functions in terms of which the existing expressions for the polyhedral 
harmonics are written. They are not, formally, the most convenient. A neater 
organisation is given at the end of this section. 
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In (22), the are the usual Legendre polynomials, in terms of which the 
usual surface spherical harmonics are, [13], 



cr(^,0) = (-ir 



{l-m)\ 



1/2 



P;"^(cos^)e*"^'^, m>0. 



(25) 



and so. 



nl/2 



2i 



{l-m)\ 
(/ + m)! 



(-1)™CJ" 



m > 
m > 



(26) 



corresponding to (24). 

The principle now to be employed, founded on (17), is that a basis for invariant 
harmonic polynomials is provided by the action of the set of independent invariant 
polynomial operators, Q{dyi, dy2, dys) = QCVy), on 1/r, (Poole, [37], Meyer, [26], 
and Laporte, [28]). This set is built algebraically from the invariant polynomial 
integrity basis. 

This technique is no different from the one mentioned above, only that, in using 
(22) the harmonic projection, (17), has really already been performed. It yields the 
polynomial series for the P[^, [35]. §85. 

In QCVy), the three-vector is effectively null. 



A = V...V,, = 2 



c>2 



= 2- 



0'' 



y y dy,dy_, dyl dy'^ dy^ dy^ 



= 0. 



(27) 



Hodgkinson and Poole derived the required integrity basis directly from the 
rotational 3-geometry of the polyhedra (the icosahedron in their case) and then im- 
posed the null-vector condition (27). Klein, [38], p. 238, calculates the polynomials 
as higher polars of binary forms (found from geometry) to which they return, up to 
a factor, on enforcing the null condition, (27), or, rather, its dual. It is, therefore, 
logically more satisfactory to start from these binary forms. ^ 

Following Honl, [36], define, in a symbolic way, the binary pseudo-operators. 



Ai 



/2 



d 



dy^ 



1/2 



d 



dy-i 



1/2 



A-1 



/2 



d 



dy-^ 



1/2 



d 
dyi 



1/2 



(28) 



^ A detailed account of the construction of the integrity bases using 3-vectors is given in Appendix 
C of [39]. 



10 



so that, in accord with (27), 



d 

^ = ±v^Ai/2A_i/2. (29) 



With these correspondances, and the upper sign in (29),^ it can easily be 
checked that ^.nd A_i/2 transform as standard spin-1/2 spinors under SU(2) 
(c/ the interesting paper by Honl, [36]) and hence invariant polynomial operators 
can be obtained from Klein's invariant binary forms on making the replacements 
Zi \ii2 and Z2 ^-1/2 fc>r his binary variables. For notational, and comparison, 
convenience, I will use z-^^ and Z2 for the operators, as well. 

Klein's fundamental icosahedral binary invariant ground form is the 12-ic 
(spin-6), 

/ = ^1^2(^1'° -^2°) +11^1 4 (30) 

in canonical form. 

This has the associated Hessian (spin-10), 

H = -(^20 + ^20) ^ 228 iz,Z2f{zl' - 4°) - 494zi°4° . (31) 

and the Jacobian transvectant (spin-15), 

T = {zf + zf) + 522 iz,Z2)'{z',' - zl') - 10005 (z.z^riz',' + zf) . 

The tesseral harmonics (24) can be rewritten, with the operator replacements, 

Niil - m)l y,^^ r'+\z,Z2y-' (zj-^ + (-l)"^zf'") - 

> (32) 
iN^il -m)\y,Y^ = r^+\z,Z2y-^ [zl-^ - (-l)-^^j - , 

where A*"; = i^2^/^/2 is an irrelevant i.e. overall, constant. 

As examples, one has the invariant harmonic forms r^/l/r, r^^Hl/r and 



^ In this case, for consistency because of (23), it is necessary to make the negative sign correspond 
to the complex conjugate, (d/dyo)* = — Ai/2^-i/2- 
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r ° T 1/r, given by, 

^ 3^c S + 228 5! % + 247 10! y, % 

~ ° - V429/13 (I?i°5° - 0) + ^46189/247 {Vl\^ + (33) 
^16 rp- r^y.ll- 522 5! II - 1000510 10! ^5 



~ (Pis^ + P^5^ °) - 87^7590/1000510 (l^i^o^ - ^^io°) 
- 3\/3338335/1000510 (©^^5° + P\^5°) , 

where the ~ sign indicates equahty up to an (irrelevant) overall constant factor. 

The expressions (33) agree with those listed by Zheng and Doerschuk, [40], 
who obtained them by a complicated method from first principles using projection. 
Making use of Klein's work in deriving his forms, many evaluations can be done by 
hand. 

The modes, (33), pull back to vertical modes on S^/F'. The remaining modes 
are obtained, in accordance with previous remarks, simply by replacing the right 
index on the Ps by m, with — Z < m < Z, for each appropriate 1. 

A somewhat similar technique is given by Kramer, [41], although he starts 
with the highest weight modes on S^, V^J (in my conventions), and lowers the / 
index. This means he does not use the puUback description and works only on the 
three-sphere. 

The general method for the construction of all polyhedral harmonics is given 
by Hodgkinson; see also Laporte [28] and Meyer, [26]. The number of independent 
invariant harmonics is the number of terms of the form in the general 

polynomial, homogeneous of degree 21 as and Z2 are equally scaled. However, 
because of the necessary syzygy, 

= 1728/^ - , (34) 

c is either zero (even, Neumann modes) or one (odd, Dirichlet modes), [28]. The 
number of modes, di, is encoded in the icosahedral generating function, [28], [26], 
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the numerator, a^^, corresponding to the odd modes ^. 

In order to expand the harmonic in any particular case, set 



T'+^rH^T--^Y.^l-m)\A'ily,^ (35) 



m=0 



with the upper hmit = 6a + 106 + 15c, c = 0, 1. It can be seen that the problem 
reduces to the algebraic one of organising the polynomial into a sum of 

the expressions appearing on the right-hand side of (32), [29]. 

For this purpose, it is suflBcient to set = 1, so that, from (32), we seek 



1=1 



~ J2 {z^^"" + (-1)""+" z^-"^) , (36) 

^^1=^ m=0 



and, as indicated, I am working up to an irrelevant overall constant. 

Evaluating the left-hand side using (30) and (31), the coefficients A'^ can be 
read off yielding the expansion coefficients in (35). In more detail, 

(z^^ - 11 - z)"" (-^2° - 228 z^^ - 494 + 228 z^ - 1)'' 

X (^30 _ ^25 _ ;loo05 z^° - 10005 z^° + 522 z^ + 1^ 

6a+106+15c 

(37) 

m=0,5,10,... 

Equation (37) for c = agrees with the general form given by Poole, [43], equn.(ll). 

Machine algebra allows one to proceed as far as required. For a given I it 
is only necessary to determine the possible values of a, and 6, an easy task, and 
then extract some polynomial coefficients by commonplace routine. An example is 
I = 30 for which c = and there are two values for (a, 6), namely (5, 0) and (0, 3). 
Correspondingly, one obtains two modes from (37). I find for these, 

~ ^'-ydl - 10! 553;J° + 151120537,^^ - 20! 13090 3^,^° 

+ 25! 69585 - 30! 134761 y^ % 



^ For the purely rotational situation, there is no obligation to refer to even and odd modes. The 
numerator simply reflects the possible existence of a single factor of T, which is related to the 
syzygy. Also, the generating function is being derived here from the integrity basis. Often, the 
generating function is derived by other means {e.g. by cyclic decomposition) and information 
about the basis obtained therefrom, e.g. Sloane, [42], [30]. 
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and 

*30^^ ~ >Vc3o+5! 6843^ + 10! 157434 3^ J° + 15! 125274603^,^^ 

+ 20!774604953^c3S + 25! 130689144 3^ - 30! 33211924 3^ . 

The modes listed in [40] are linear combinations of these. Thus 

Tgo^i ~ 12251 ^'^o'^^ - 33211924/11 

7^30,0 ~ - 45750 *i'o'°^ • 

The procedure for the cubic groups is similar. The fundamental ground form 
can be taken to be the special sextic, 

f = z,z^{zt - zt) = 4 (38) 

(corresponding to the Cartesian, y^y'^y^) from which the complete form system can 
be derived, geometrically, e.g. [44], or algebraically by invariant theory, [45]. 
The Hessian, H, and the Jacobian, T, of / and H are 

H = if, ff = (abfaibi = -^{4 + 14 + 4) 
T={f,H) = -^ {zl^ - 33 zfzt - 33 zfzl + zf) , 

where (/, g) is the transvectant of the forms / and e.g. [46]. 

The form / is an absolute invariant for the tetrahedral group and so a com- 
plete set of absolute invariants is /, and T, corresponding to spins 3, 4 and 6, 
respectively. 

The standard syzygy, 

t2 = -— (39) 

108"' 2 ' ^ ' 

means that the general term can be written /" where c = 0, 1 and I = 

3a + 46 + 6c so that the generating function is. 

By contrast, the octahedral group replaces / by ±/ which means that the 
Hessian, H, is an absolute invariant, but that T is replaced by ±T. Therefore a 
basis set of absolute invariants in this case is provided by H and /T, with spins 
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6, 4 and 9, respectively, e.g. [38], p. 69. Comparing notations with Laporte, [28], 
f ^ (f)^, H (f)^., T ~ $g, fT ~ $9, as can be checked algebraically. 
The syzygy, (39), is recast as, 

^fTf = -f(^^^f+\H'), (40) 

so that the term in the general polynomial this time is {fTY with c = 0, 1 

and Z = 6a' + 46 + 9c. The octahedral generating function is therefore, 

^(-; o) = (1 _ ^4|(f_ ,6) = E dii; o) , 

given previously. 

As a typical example set / = 12 so that c = and (a', 6) = (2, 0) or (0, 3). The 
general definitions (35) and (36) still apply with the upper limits = 3a + 46 + 9c, 
where a = 2a' and, in place of (37), one has 



(^ - z^Y (1 + 14^4 + z^YiX - 33^^ - 33z« + z^'^y 

3a+46+9c 

= J] (^^+"^ + (-1)"^+= 2'-"^) , 



(41) 



which readily yields the invariant modes, (c = 0), 

*Sf^~4!3;,«,-8!3;,^2 + 12!33^ei2 

and 

~ :^^c + 4! 42 ?2 + 8! 591 \^ + 12! 1414 \^ . 

These quantities agree with the corresponding ones in Table 8 in Altmann 
and Bradley, [47], after a linear combination which orthogonalises them. I have 
not checked them against the Cartesian expressions given by Ronveaux and Saint- 
Aubin, [32], equn.(34). 

Incidentally, it is interesting to note that the fourth transvectant of the octahe- 
dral form / with itself is zero. In fact, this characterises the form, (38), Klein, [44], 
Cayley, [48], Gordan, [45]. In angular momentum terms this means that the van- 
ishing of the Clebsch-Gordan spin-two combination of two spin-three quantities, 

\ I ran I 
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implies that / is equivalent to the special octahedral form, (38). 

Similar remarks hold for the icosahedral case and are related to Fuch's notion 
of Primformen. 

A further ancillary point concerns the nature of the syzygies. These can be 
proved either geometrically, [38], or algebraically, [45]. Another approach uses the 
action under the reflections of the extended groups. Taking the icosahedral syzygy, 
(34) as an example, T is odd under reflection while H and / are even. Hence T^, 
being even, can be expressed algebraically in terms of / and H, e.g. [28]. In fact, all 
the left-hand sides of the syzygies, (34), (39) and (40) are squares of odd quantities 
(Jacobians) and the corresponding modes thus are counted by the numerator in the 
generating functions given above. This is a well known behaviour. The denomi- 
nator basis elements sometimes are referred to as 'free', and the numerator ones 
as 'constrained' and occur only once in the construction of the algebraic basis, e.g. 
Cummins and Patera, [49] , Patera et al, [30] , McLellan, [50] . 

Of course the entire scheme is a classic example of the invariants of finite 
reflection groups and can be treated from this point of view ah initio. 

For the cyclic and dihedral groups, the construction of invariant bases in terms 
of the Legendre functions proceeds simply by a process of selection and is given in, 
e.g.., Meyer, [26] and, earlier, in the classic, Pockels, [51]. See also Altmann and 
Bradley, [47]. 

In order to make comparison with existing results easier, I have used the tesseral 
harmonics. However it is more elegant to use the usual spherical harmonics, (25). 
Then, instead of (32), some rearrangement yields the neater relation, 

KO 1 
r 

with the familiar monomial (or null Z-spinor), 



{-iyN,CUeA)=r'^'^^^-. (42) 



/ + n\ I — m 



{{I - m)lil + m)\) 



1/2 



constructed from the two-spinor ( ^ ) , c/ [36], equn.(28). 



However, we have now come full circle, because (42) can be rewritten 

||^C^(^,0) = r'+^Ct,(-V,)i, (43) 



where C^(a) is the solid spherical harmonic and is a null vector according to 
(27). Equation (43) is an example of Nivens' general theorem, [52], [35], p. 127, and 
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could be taken conveniently as the starting point of the analysis, rather than the 
specific (22). It also follows from (17). 



6. Tensor polyhedral modes 

The modes found in the last section are invariant, i.e. they transform as the 
identity rep, 1, of T. In order to find those that transform equivariantly according 
to the other irreps (some are listed in [47] for example) I use the results of Patera 
et al, [30], and, especially, of Desmier and Sharp, [53] and of Bellon [8]. I restrict 
myself to the octahedral and icosahedral groups. 

To start with, the mode (which could be called a twisted scalar mode) trans- 
forming as the one-dimensional rep, l', of the octahedral group has the general 
form 

where A and B are invariants constructed as in the preceding section from and 
H. For example, for Z = 12, ^ = and B p and the mode is easily calculated 
to be, 

^12(1') ~ 10! Uy^^^ - 6! 35 + 2! 37,}° , 

again in agreement with [47] Tabic 8. 

The next most complicated modes transform according to the irrep 2 and have 
the general structure 

*,(2)=r'+i^,(2)i 

where, [53] (18), 
^^(2) = 

VZ{zt + Qzlzl + z|) \ / V^{zl + Azlzl - lOzfzl + Azfzl + zl) 

with A and B as before. Choosing again I = 12, and writing 

A^f'H^ B^P^'h^', 

the only possibilities are (a, h) = (1, 1) and (a', b') = (0, 2) giving two independent 
modes. I write them out just to show how the method proceeds. 



17 



Applying the technique explained in the last section directly to the 2-tensor 
integrity basis given by Desmier and Sharp yields the two modes, (I write y for 
for space reasons), 



V3(-12!783^02 - 10! 143^^2 + 81723^4^ + 6!l33^f2 + 4!63^f2 + '^^^12. 
-12!783^?2 + 101423^2^ + 81723^4^ - 6!393^f2 + 4!63'f2 - '^^■^yl2 



(44) 



and 



/\/3(- 1219623^02 + 1019043^2^ - 8!813^i4^ + 611163^^2 + 41183^^2 + 2!3^i^2°) 
*?2(2) = ( - 1219623^02 + 10!27123^i2^ - 81813^^2 - 6!3483^f2 + 41183^^2 

V -m2ylUyE) 

(45) 

which are, in appearance, more complicated than those in [47] Table 10. The reason 
is that the basis in 2-irrep space used in [53] is different to that in [47]. 

In order to convert the modes (44) and (45) to Altmann and Bradley's, firstly 
multiply them by the S0(2) rotation. 



R 



cos7r/3 — sin7r/3 
sin7r/3 cos7r/3 



to convert the 2 basis, and then take linear combinations in order to get precise 
agreement. I find 

■^^i\2) ~ 37i?*i2(2) - 3i?*?2(2) 
'(2)~8i?*}2(2) + i?*?2(2). 

The same considerations apply to the other irreps of the octahedral (and the 
binary octahedral) group but I do not need to give any further examples. I comment 
that the 3-irrep of O (r4 in [53] ^O) requires no rotation of basis. For example the 
fundamental I = 1 mode is, 

f(2) ^ 

with the notation of [53], and so has a Cartesian basis. This agrees as it stands 
with the entry in [47], Table 11. Tg of [53] is the spinor, quaternion irrcp, 2^.. 

The general structure of the bases is outlined in Section 4 of [53] and it is 
clear that any desired mode can be constructed without too much bother. A basic 
symbolic manipulator (I used DERIVE) is all that is required. If one were to pursue 




There appears to be a misprint in equn.(56) of [30]. The irrep for O should be 
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this extensively, it would be advisable to rotate the 2-octahedral expressions given 
in [53] by R at the start. 

The generating functions for the icosahedral group are given in Table VII of 

[53] but not the corresponding reps. However, Bellon, [8], gives a list of those based 
on the fundamental quaternion irrep, 2^, which can be used in exactly the same 
way as above to give the corresponding tensor icosahedral harmonics. 

A technical point is that the reps in [8] are taken with respect to a spherical 
basis and it is easiest to use standard spherical functions and the basic relation (42). 

As an illustration, the I = 1 and / = 5 modes for the vector 3-irrep are rapidly 
found to be proportional to. 



respectively. In this form, the lifting to dodecahedron space is easy. 

It would need a certain motivation to compute further cases or to check or- 
thogonality for degenerate modes, although there is no particular difficulty. 

The modes which transform according to the spinor irreps of the binary group, 
r', will involve what might be termed polyhedral spinor harmonics leading to poly- 
hedral spinor hyper harmonics on S^/V. These will be dealt with at another time 
using fractional derivatives. 

7. Concluding remarks 

I have employed the Hopf map to lift modes on the two-sphere to modes on 
the three-sphere, a known procedure, and have divided by a polyhedral symmetry. 
I used this to motivate the construction of the required symmetry adapted modes 
on the orbifolded two-sphere which are computed by a binary method and agree 
with those obtained many years ago by the more standard, and in my view more 
involved, technique of projection. This has been extended to tensor modes under 
the action of the deck group and some specific cases have been evaluated using the 
results of Desmier and Sharp and of Bellon. 

Expressing the symmetry adapted modes in the traditional way as sums of 
spherical harmonics allows the complete set of modes on the factored three-sphere 
to be found from the puUbacks using raising and lowering with no extra work. 




v 
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An aspect that deserves exploration is the significance of the degeneracies, 
d{l; r), as the dimensions of the irreps of the symmetry group of S^/F which is the 
centrahser of F in S0(3), or for S^/F', the centrahser of the binary group, F', in 
SU(2). This is the 'missing label' question. 
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